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In the presence of a uniform field the one-dimensional spin-± antiferromagnetic Heisenberg model 
develops zero frequency excitations at field-dependent 'soft mode' momenta. We determine three 
types of critical quantities, which we extract from the finite-size dependence of the lowest excitation 
energies, the singularities in the static structure factors and the infrared singularities in the dynam- 
ical structure factors at the soft mode momenta. We also compare our results with the predictions 
of conformal field theory. 



PACS number: 75.10 -b 



> 

in 
oo 

o 
^t- 
o 

Q\ ; I. INTRODUCTION 



In this paper we are going to study the zero temperature dynamics of the one dimensional spin-^ antiferromagnetic 
Heisenberg model 



N N 



H = S(x)S(x + 1)-2BJ2 S 3 (x), (1.1) 

x — 1 x—1 

in the presence of a uniform external field B. The quantities of interest are the dynamical structure factors at fixed 
magnetization M = S/N: 
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S a (u;,p,M,N)=Y / 5(u-{E n ~E s ))\(n\S a (p)\s)\ 2 , a = 3,+,-. (1.2) 

n 

They are defined by the transition probabilities | (n\S a (p)\s) | 2 from the ground states |s) = \S, S3 = S) in the subspaces 
with total spin S and energy E s to the excited states \n) with energy E n . The transition operators we are concerned 
with, are the Fourier transforms of the single-site spin operators S a (x): 

1 N 

S a (p) = —J2e ipx S a (x), a = 3,+,-. (1.3) 

V X— 1 

The structure factor s fll.2| ) have been investigated before by Muller et alM They performed a complete diagonalization 
of the Hamiltonian ( |1.1| ) on small systems (N < 10) and analysed the spinwave continua by approximately solving 
the Bcthe ansatz equations for the low lying excitations. In particular, they found a lower bound 
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uo> \u 3 (p,M)\, 



(1.4) 



W3 ( P ,M) = 2i ? sin|si n ^W ; 



(1.5) 



for the exci tations contributing to the longitudinal structure factor Ss(u!,p, M). The constant D on the right hand 
side of (1.5) is fixed by the magnetization curven 



B(M) = 2D$,unrM. 
The lower bound vanishes at p = and at the field dependent momentum 

p 3 (M) = tt(1 - 2M), 



(1.6) 



(1.7) 



signalling the emergence of zero frequency modes ('soft modes') in the spectrum of excitation energies. The analysis 
of the spinwave continua relevant for the transverse structure factors S±(u),p, M) leads to the following approximate 
lower bounds 



w > ^±(p, M), 

for the excitations produced by the raising and lowering operators S+(p), S-(j>), respectively: 



lu + (p, M) = 2D sin | cos (| - ttM) - simrM 



for pi(M) < p < it, 



and 



w-(p,M) = \w 3 (ir-p,M)\ for < p < it. 



(1.8) 



(1.9) 



(1.10) 



Both bounds vanish at p — it and at p — p\(M) — 2ttM. The softmodes at the field dependent momenta pj(M), jL= 
1,3 produce characteristic structures in the momentum dependence of the corresponding static structure factors E'u 
It is the purpose of this paper to an alys e the singularities in the static structure factors and the infrared singularities 
in the dynamical structure factors (1.2) at the softmode momenta. In Sec. II we review our method to compute the 
excitation energies and transition probabilities for finite rings (N < 36). The finite-size dependence of the lowest 
excitation energy at the soft mode momenta is analysed by solving the Bethe ansatz equations on large systems 
(N < 2048). The critical behavior of the static structure factors at the softmode momenta p = p a (M), a = 1, 3 and 
fixed magnetization M — 1/4 is investigated in Sec. Ill based on a numerical computation of the ground state on 
rings with N = 12, 16, ...,32,36 sites. In Sec. IV, we demonstrate how infrared singularities emerge in a finite-size 
scaling analysis of the dynamical structure factors in the euclidean time representation. Finally, in Sec. V we compare 
our numerical results with the predictions of conformal field theory. 



II. SOFTMODES IN THE EXCITATION SPECTRUM. 



An approximate scheme to determine low lying excitation energies and transition probabilities has been proposed in 
Ref. [ 0. It starts from the recursion algorithm,B which generates a tridiagonal matrix. Eigenvalues and eigenvectors 
of this matrix yield the exact excitation energies and transition probabilities. There are, however, two sources for 
numerical errors in this scheme. The orthogonality of the states produced by the recursion algorithm is lost more and 
more with an increasing number of steps, due to rounding errors. Moreover, the iteration has to be truncated before 
the Hilbert space is exhausted. 

Nevertheless the method yields good results for the lowest 10 excitations - provided that these contain the dominant 
part of the spectral distribution. This condition is satisfied for the excitations in S a (u>,p, M, N), a = 3,+. For 
S~(uj,p, M, N) near the softmode momentum p\(M), however, this is not the case. In Table I we compare the low 
lying excitations for S-(u),p, M, N), M = 1/4, p = tt and p = tt/2 — 2ir/lQ on a ring with N = 16 sites, as they follow 
from an exact diagonalization (upper part of Table I) and the recursion algorithm (lower part of Table I) , respectively. 

At p = 7T, 76.95% of the spectral weight is found in the first excitation. Energy and relative spectral weight of 
the first excitation are reproduced within 13 digits. The following 7 excitations can be identified term by term with 
decreasing accuracy for the energies and the relative spectral weights. 

The situation is different forp_ = tt/2 — 27r/16, which can be seen in the right hand part of Table I. The exact result 
yields large spectral weights - marked by an asterisk - for the 1st (19.55%), the 15th (18.33%) and the 20th (13.80%) 
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excitation. The recursion method reproduces the energy and spectral weight of the first excitation within 13 digits. 
The two other excitations with large spectral weight - marked by an asterisk - are only in rough agreement with the 
exact result. We fou nd, h owever, that this inaccuracy has no effect on the dynamical structure factors in the euclidean 



time representation (4.1). The latter will be investigated in Sec. IV. In Figs. l(a),(b), (c) we present the momentum 
dependence of the excitation energies in the dynamical structure factors S a (ui,p,M — 1/4, N — 28) as they follow 
from the recursion method. The size of the symbols measures the relative spectral weight \ (n\S a (p)\s}\ 2 / S a (p, M, N). 
The normalization is given by the static structure factors: 

/>oo 

S a (p,M,N)= / du>S a (u,p,M,N), a = 3,+,-. (2.1) 

Ju a (p,M,N) 

There is a strict relation between the static transverse structure factors: 

S-(p, M, N) = S + (p, M, N) + 2M. (2.2) 

It should be noted that S+(p, M, N) ss for p < pi(M) [cf. Fig. 3(b)], which implies that the absolute spectral 
weight |(n|5+(p)|s)| 2 is almost zero for p < pi(M). 



The solid curves represent the lower bounds ( |1.5| ),( |l.9| ) and ( 1.1C ) obtained from the analysis of the spinwave 
continua.El The emergence of the softmode at p — p 3 (M = 1/4) == n/2 in the longitudinal c ase [Fig. 1(a)] is clearly 
visible. Note, that there are some excitations with small spectral weights below the bound ( |l.5| ) (for p > 3%/ 4). We 
do not know, whether the spectral weights will survive in the thermodynamical limit. 

The lowest excitations in the transverse cases [Figs. 1(b) and 1(c)] are found at p = ir and at the field dependent 
momenta 

2ir 

pt(M)= Pl (M)±—. (2.3) 

We have analysed the finite-size dependence of the lowest excitation energies: 

u 3 (p 3 (M), M, N) = E(p = Ps + p 3 (M),M = S/N, N) - E(p s ,M = S/N, N), (2.4a) 
- ,(-. M. \ ) = E(p = p s + n, M = (S + 1)/N,N) - E{p s ,M = S/N,N), (2.4b) 
lo ± (p = pf(M),M, N) = E(p s +pf(M),M ={S± 1)/N, N) - E(p s ,M = S/N, N), (2.4c) 

p s denotes the groundstate momentum in the sector with total spin S; p s = if N + 2S is a multiple of 4, p s = 
7r otherwise. The lowest energy eigenvalues E(p, M, N) with momentum p and spin S were computed on large 



systems (N < 2048) by solving the Bethe ansatz equations. The extrapolation of the energy differences (2.4) to the 
thermodynamical limit 

Urn Nlj 3 (p 3 (M),M,N) =fi 3 (M), lira Nui(n, M,N) = fii(M), (2.5a) 

Um Nw± (pf (M) ,M,N) =flf(M), (2.5b) 

obey the following relations: 

Qf(M) = n 3 (M) ± Qi(M). (2.6) 

Together with the spinwave velocity v(M) 

2nv(M) = lira N[E{p s + 2tt/N, M, N) - E(p s ,M, N)], (2.7) 

they define the scaled energy gaps: 

29 a (M) = ^^±, a = 3,1, (2.8a) 

2flf(M) = = 2[0 3 (M) ± 6\{M)]. (2.8b) 
The M-dependence of the quantities 8 a (M), a — 3,1, is shown in Fig. 2. It turns out that 
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2 * l(M) = 2UM) ^ 

in accord-jwith the analytical result of Bogoliubov, Izergin and Korepinl. In the limit M — > 1/2 one finds 2# 3 (M) = 
1 + 2MJB The dotted line in Fig. 2 near M = indicates the logarithmic singularity 



20 3 (Af)^l+(Jn— j , (2.10) 
which was obtained by Bogoliubov, Izergin and Korcpin0 by a perturbative approach to the Bethe ansatz equations. 



III. CRITICAL BEHAVIOR OF THE STATIC STRUCTURE FACTORS AT THE SOFTMODE 

MOMENTA 

The static structure factors of the antiferromagnetic Heisenberg model in the presence of a magnetic field have been 
investigated in a previous numerical study on systems up to N — 28 u Meanwhile we have extended the system size 
to N = 32 and N = 36 at fixed magnetization M = 1/4. We find the following features: 

1. The transverse structure factor at momentum p = n diverges for N — > oo. A power law fit 

Sx(Tt,M,N) 0.5037V 1 -" 1 (M) , (3.1) 

to the finite system results for N = 36, 32, 28 leads to the value rji(M = 1/4) = 0.65 for the critical exponent. 
The same exponent governs the approach to the singularity in the momentum p: 



p — >TT _ . _ P 



7T 



Vl(M) — l 



Si(p,M,oo)^ 0.316 1-- , (3.2) 



The finite-size dependence (3.1) is shown in Fig. 3(a). The momentum dependence can be seen in Fig. 3(b) 
where we have plotted S\(p = tt,M = j,N) versus (1 — £) T >i( M )~ 1 using the critical exponent determined in 
Fig. 3(a). 

2. The approach to the field dependent soft mode pi(M) = 2ttM in the transverse structure factor is shown in the 
upper left [p — > pi(M) — 0] and the lower right [p — > p±(M) + 0] insets of Fig. 3(b). The numerical data behave 
as 



Si(p^pi(M)±0,M,oo) 



P 



Pi(M) 



nf(M)-i 



(3.3) 



if the critical exponents are chosen to be T)i(M = 1/4) = 2.17 , %(M — 1/4) = 0.8. ..1.2. The uncertainty 
i n rf x(M = 1/4) reflects an instability in the fit to the numerical data. Note, that the right hand side of 
( |3.3[ ) diverges for T]i(M — 1/4) < 1 but converges for T]i(M = 1/4) > 1. An unambiguous determination of 
rf[{M = 1/4) demands for much larger systems than N = 36. 

3. The finite-size dependence of the longitudinal structure factors at p — p${M) 

S 3 (p 3 (M), M, N) N -^f -Q.l2AN l -^ M ) + 0.308, (3.4) 

is shown in Fig. 4(a) for M — 1/4, p = p$(M) = ir/2. A power law fit to the finite system results with 
N = 36,32,28 yields: r] 3 (M — 1/4) = 1.51. The same exponent governs the approach to the singularity from 
the left: 

/ \»?3(Af)-l 

S 3 (p^p 3 (M)-0,M,N) N ^-0m2\l--^—) +0.322, (3.5) 

as is demonstrated in Fig. 4(b). It is not so easy to decide, whether a different exponent is needed to describe 
the approach to the singularity from the right. In the inset of Fig. 4(b) we plotted the approach from the right 
versus |1 — p/p 3 (M)\ r i 3 ^ M=1 ^~ 1 . 
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The Fourier transform of the singularities in the static structure factors determines the large distance behavior 
of the corresponding spin spin correlators: 



(s\S 1 (0)S 1 (x)\s) ™ 
(s\S 3 (0)S a (x)\a) - (s\S 3 (0)\s) 2 co S [p 3 (M)x 



. . A\ (M) , , . . / Af(M) A7(M) 
cos(7ra) 1 \' + cos[ Pl (M)x' 1 1 



x m(M) 

MM) 
x m(M) ■ 



(3.6a) 
(3.6b) 



Confo rma Uicld theory! predicts a relation between the critical exponents r)(M) in (3.6) and the scaled energy 
gaps j2^):BB 



26 a {M) = j] a (M), a = 3,1, 
20f(M) = t]f(M). 



(3.7a) 
(3.7b) 



A derivation of (3/7) is presented in appendix A. A comparison of the left and right hand sides of (3/7) is 
presented in Table II. 



IV. FINITE-SIZE SCALING ANALYSIS OF THE INFRARED SINGULARITIES 

The euclidean time representation 
S a (r,p,M,N) -- 



dtu e-" T S a (w,;p,M,A0, a = 3,+,-, 



(4.1) 



ui a {p,M,N) 



is most suited to study finite-size effects in the dynamical structure factors ([L^). The singularities in the static 
structure factors S a (r = 0,p, M, N) at the softmode momenta originate from the infrared singularities in the dynamical 
structure factors. In the combined limit 



N -> oo, 



keeping fixed the 'scaling' variables 



z a (p,M) = TUj a ( P ,M,N), a = 3,+,-, 



(4.2) 



(4.3) 



the low frequency part at the softmode momenta p — ir,p = pi(M) ± 2ir/N,p = p 3 (M) is projected out. We therefore 
expect to see here directly signatures for the infrared singularities. Let us assume that the emergence of the infrared 
singularities on finite systems can be described by a finite-size scaling ansatz: 

S a (u,p,M,N)=u- 2a «^g a (u;/u) a ( P ,M,N),n a ( P ,M,N)), a = 3,+,-. (4.4) 

The scaling functions g a are supposed to depend only on the scaled excitation energies u/L0 a (p, M) and the variable 

n a (p, M, N) = [p- p a (M)]N/(2ir), (4.5) 

which describes the approach to the softmode m omenta. The ansatz ( kl.4 | ) i nduces the following finite-size scaling 
behavior of the euclidean time representation ( |4.1| ) in the combined limit ( pl.^ ) and (4.3): 



r i-2aab,M) 5a(r)P) Mj N) = Ga ( Za ( P: M),n a (p, M, N)) cxp[-z a (p, M)]. 
The two scaling functions on the right hand sides of equations (fh4|) and (fb6j) are related via: 

i 

dxe- {x - 1)z g(x,n). 



G(z, n) 



A-2a 



N = 16, 20,. ..,36 and a: 



(4.6) 



(4.7) 



N = 16, 20,..., 32 



Based on our numerical results for S a (r,p, M, N) at M = 4, a = i 
at the softmode momenta we will now test the validity of the finite-size scaling ansatz ( |4.6| ). 

Let us start with the longitudinal structure factor at the s oftm ode p = p${M = 1/4) = n/2. In this case the variable 



(4.5) is n 3 (p = n/2,M = 1/4) = 0. The left hand side of (4.6) versus the scaling variable z 3 (p = w/2,M = 1/4) is 
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shown in Fig. 5(a) for the following values of a 3 (p = n/2,M = 1/4) = 0.22,0.23,0.234. For z 3 > 0.4 (inset of Fig. 
5(a)) the finite system results coincide best if 

a 3 (p = tt/2, M = 1/4) = 0.23. (4.8) 

Therefore, this is the expected critical exponent for the infrared singularity in the longitudinal structure factor. 



Deviations from this value for a 3 on the left hand side of (4.6) obviously lead to a violation of finite size scaling. It 



is remarkable to note that finite-size scaling [with the exponent a 3 (p = ir/2,M = 1/4) = 0.23] persists for all values 
z 3 > 0.4. In the limit z 3 — > oo the first excitation alone survives and we can conclude for the finite-size dependence 
of the transition probability: 

\(n = l\S 3 (p = vr/2)| S )| 2 ^-=5° N 2 ^- 1 . (4.9) 

In other words, the critical exponent a 3 for the infrared singularity can by read off the finite-size depeadence of the 
transition probability for the first excitation. Indeed this feature is predicted by conformal field theory.Q (cf. (A9) in 
appendix A) 

Next we turn to the infrared singularities of the transverse structure factors S±(u>,p = tt, M — 1/4). As can be 
seen from Fig. 5(b), finite-size scaling is found for the following choice of the critical exponents: 

a + (p = n,M = 1/4) = 0.69, (4.10a) 
a-(p = n,M = 1/4) = 0.66. (4.10b) 

In contrast to the longitudinal case, finite-size scaling can be observed here for all values of the scaling variables 
z + ,z_. 

Finally we present in Figs. 6(a) and 6(b) the tests of finite-size scaling for the transverse structure factors S±(r,p = 
tt/2 ± 2it/N, M = 1/4, N) if we approach the field dependent soft mode p x {M = 1/4) = tt/2 from the left (p = 
tt/2 — 2tt/N) and from the right (p = tt/2 + 2ir/N), respectively. The critical exponents are found to be 

a+(p = tt/2 + 2n/N, M = 1/4) = -0.20, (4.11a) 
a_(p = tt/2 - 2n/N,M = 1/4) = -0.05. (4.11b) 

Finite-size scaling works quite well for S+ for large and small values of the scaling variable z+ as can be seen from 
the inset in Fig. 6(a). This is not the case for S-. Here finite-size scaling breaks down for small values of z_ as is 
demonstrated in the inset of Fig. 6(b). The critical exponent a_ (p = tt/2 — 2tt/N, M = 1/4) = —0.05 results from the 
finite size scaling analysis for large values of z_, where the transition probability for the first excitation is projected 
out and has the following finite size dependence: 

|(n= l|S*_(p = 7r/2-27r/A0|s)| 2 N 2a — 1 . (4.12) 



V. DISCUSSION AND CONCLUSIONS. 



In the presence of a uniform field, the one-dimensional antiferromagnetic Heisenberg model is critical in the following 
sense: The excitation spectrum is gapless at the momenta p = 0, p = tt, p = p 3 (M) = n(l — 2A1) and p — pi(M) = 
tt ■ 2M . In this paper we have tried to answer the following question: Is conformal field theory applicable to describe 
the low energy excitations at these momenta ? To answer this question we have determined: 



1. the scaled energy gaps 26(M), defined through (|2.4|)- ( 



2. the critical exponents rj(M) for the singularities (3.2), (p.3|) and (3.5) in the static structure factors 



3. the exponents a(M) for the infrared singularities (4.4) in the dynamical structure factors 



A compilation of the various critical quantities for M = 1/4 is given in Table II. 

The predictions of conformal field theory are reviewed in appendix A. In particular the following relation is expected 
to hold: 



26{M) = jj(M) = 2[1 - a(p,M)} 



(5.1) 



Looking at Table II we find: 
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(a) The critical quantities 28 3 (M = 1/4), i] 3 (M = 1/4) and 2 - 2a 3 (p = tt/2, M = 1/4) agree within the numerical 
uncertainty. Moreover, the critical exponent a 3 (p — tt/2, M = 1/4) also governs the finite-size dependence of 
the transition probability for the lowest excitation (|0^). We therefore conclude, that the excitations in the 
longitudinal structure factors at the softmode p^(M) = 7r(l — 2M) are correctly described by conformal field 
theory. 

(b) The critical quantities 2&i(M = 1/4), r)i(M = 1/4), 2 - 2a + (p = tt, M = 1/4), 2 - 2a-(p = tt,M = 1/4) agree 
within numerical uncertainties. In both cases the finite-size dependence of the transition probability for the 
lowest excitation is in accord with the prediction of conformal field theory. 

(c) The critical quantities 28f(M = 1/4) and r]f(M = 1/4) agree within numerical uncertainties and deviate by 
about 15% from the exponent 2[1 - a + (p = tt/2 + 2n/N,M = 1/4)]. 

(d) The scaled energy gap 28 ^ (M — 1 /4) agrees with the critical exponent rf[ (M = 1/4) - within the large numerical 
uncertainty - but strongly deviates by more than a factor of 2 from the exponent 2[1 — a_ (tt/2 — 1/(2 AT), M — 
1/4)], which we extracted from the finite-size scaling analysis of the infrared singularity in the transverse 
structure factor S- at the softmode p = pi(M) — 2tt/N, M = 1/4. It was demonstrated in Fig. 6(b) that finite- 
size scaling only works for large values of the variable z_, where the first excitation alone contributes. Therefore, 



the exponent 2[1 — a(jr/2 — 2tt/N, M — 1/4)] is fixed by the finite-size behavior (4.12) of the transition probability 



for the first excitation. The exponent is definitely different from the scaled energy gap 28 \ (M — 1/4). 
It is worth while to note that in the cases (a), (b) and (c), where we find agreement of our numerical results with the 



prediction (5.1) of conformal field theory the spectral weight of the excitations is concentrated at low frequencies. 
This can be seen directly for the case (b) (p — tt) in the left hand part of Table I. In contrast, the right hand part of 
Table I shows the widespread distribution of the spectral weight for case (d). Here we were not able to establish the 



identity (5.1) 
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APPENDIX A: CRITICAL EXPONENTS IN CONFORMAL FIELD THEORY 

In the absence of a magnetic field the spin-i Heisenberg model is known to be conformal invariant. Switching on 
the magnetic field rotational invariance is broken explicitly. Nevertheless the system remains gapless. Let us assume 
that the low energy physics of the model is still governed by conformal field theory. Then the dominant contribution 
to the long distance asjjmptotics of the zero-temperature dynamical correlation functions in the infinite x — t plane is 



correctly described as£2 

( S | S „(0,0 )S „(M)W - W4(M) W » - ^■ <M \ x + v{M}t? ^. v{M}t? K^r < A1 > 



v(M) is the spin wave velocity defined in (2.7), A a (M) and A a (M) are the conformal dimensions of the operator 
S a (x, t). The dynamical structure factor S a {u,p) is just the Fourier transform of flAip with an appropriate regulariza- 
tion. The latter can be achieved by giving an infinitesimal imaginary part to the spinwave velocity v(M). Standard 
methods yield 

S a (u,p) ~ {uJTv{M)[p-p a {M)]} 2 ^ M)+2A « {M) -\ (A2) 

near the singularities 

to ^ ±v(M)[p - Pa (M)]. (A3) 



Equation (A2) is obtained if we first consider the case A a (M) + A a (M) > 1/2 and then continue analytically. A 



conformal transformation to a strip geometry of width A^ tells us how the conformal dimensions A a (M) and A a (M) 
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are related to the energy and momentum of the lowest excitation |1) provided that the transition matrix element 
(s|5 a (0, 0) 1 1) does not vanish: 



2A a (M) = 9 a (M) + n a , 
2A a (M) = a {M) — n a , 



where 



la = [p-Pa{M)} 



N 
2tt' 



Therefore we conclude that the infrared singularity of the dynamical structure factor 

1 



is independent of n a : 



{uo±v{M)[p-p a (M)]} 2a « {My 
a a (M) = l-e a {M). 



The critical exponent rj a (M) can be read off directly from (Al): 

Va (M) = 2A a (M) + 2A a (M) = 29 a (M). 



(A4a) 
(A4b) 



(A5) 



(A6) 



(A7) 



(A8) 



In (Al) it is assumed that the coefficient A a (M) is nonvanishing. From the conformal transformation to the strip 
geometry a relation between A a (M) and the transition matrix element can be derived: 



A a (M)= lim 



N 



28 a. (M) 



<\(s\S a (x,0)\W 



Therefore, the matrix element is expected to scale as 

\(s\S a (x,0)\l)\ 2 ~ N 2a «W- 2 . 
If a finite-size analysis of these critical exponents reveals that 

6 a (M) <l-a«(M), 



(A9) 



(A10) 



(All) 



the coefficient A a {M) vanishes. In this case the expression (Al) does not represent the dominant contribution to the 
dynamical structure factor. 
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Figure Captions 

FIG. 1. Momentum dependence of the excitation energies in the dynamical structure factors at M = 1/4: (a) 
S 3 {u),p,M = 1/4, N = 28), (b) S+(uj,p,M = 1/4, AT = 28), (c) 5_(w,p,M = 1/4, AT = 28). The relative spectral weight 
is characterized by the different symbols. 



FIG. 2. The dependence of the scaled energy gaps 9i(M), 9s(M) on the magnetization M. 

FIG. 3. The transverse static structure factor at M — 1/4: (a) finite-size behavior at p = tt; (b) the momentum dependence 

(I-p/tt)" 1 ^) -1 forp^yr; (1 - 2p/7r)"r (M)_1 for p -> tt/2 - (inset upper left) and |1 - 2p/7r|^ (M)_1 forp^7r/2 + (inset 
lower right) 



FIG. 4. The longitudinal static structure factor at M = 1/4; (a) finite-size behavior at p = pz(M) = tt/2 (b) the momentum 
dependence: j 1 — 2p/7r I'wW-i for p < n/2 and p > 7r/2 (inset), respectively. 

FIG. 5. Test of finite-size scaling for the infrared singularities in the dynamical structure factors at M — 1/4: (a) the 
longitudinal case at the softmode p = ps(M) — tt/2. The inset resolves scaling violations for small values of the scaling variable 
z+; (b) the transverse cases at p = tt. 

FIG. 6. Test of finite-size scaling for the infrared singularities in the transverse structure factors at M = 1/4: (a) the 
transverse case S+ at the softmode p = pj(M) = tt/2 + 2n/N. The inset shows a magnification for small values of the scaling 
variable z + ; (b) the transverse case S- at the softmode p = p 1 "(M) = 7r/2 — 2tt/N. The inset resolves scaling violations for 
small values of the scaling variable Z- . 
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TABLE I. Energies and transition probabilities for the lowest excitations in the transverse structure factor S-(u>,p, M, N) 
for M — 1/4, N — 16, p = tt (left hand part) p_ = tx/2 — 27r/16 (right hand part). The upper and lower parts in the table 
contain the results of an exact diagonalization and the recursion method, respectively. 



S-(t = 0,p = tt) = 2.52360427892220 g_(r = 0,p = p_) = 5.01384876969894 ■ IP' 1 



U!„(tv) 


w„{-k) 






LO n (p-) 


tt)n(p-) 








0.24490318120407 


7.69543336339913 


10" 


-1 


0.87610327625377 


1.95498761012465 


10" 


-1 




2.00062423661784 


9.81468201828658 


io- 


-2 


2.50939624323648 


5.59400064383486 


10" 


-2 




3.16271478820513 


2.64572507440814 


10" 


-4 


3.47398478523209 


6.94575828976292 


10" 


-3 




3.57865017174411 


6.85304507352309 


10" 


-3 


3.60324922819252 


9.82162858347871 


10" 


-4 




3.98061972078759 


4.71390119166436 


10" 


-2 


3.71327071028290 


8.01415946306466 


10" 


-2 




4.35269652499191 


9.62711159680598 


10" 


-5 


4.17033985462645 


7.71455541891509 


10" 


-2 




4.72994384264668 


8.68880923919877 


10" 


-4 


4.21405414829430 


2.86564726032145 


10" 


-4 




5.11224598930047 


4.94692033004724 


10" 


-4 


4.30616024321460 


2.50148321929865 


10" 


-2 




5.25995835463119 


3.77843015439737 


10" 


-5 


4.39960077459270 


1.18063340236613 


10" 


-3 




5.45318695502460 


7.32076078229036 


10" 


-3 


4.77941756257073 


7.94747197100013 


10" 


-3 




5.74223821730404 


3.48697863610770 


10" 


-2 


4.99153366093631 


5.41969707921773 


10" 


-6 




6.14223417771473 


2.73940664197773 


10" 


-6 


5.10045741321637 


5.50730490279844 


10" 


_2 




6.20371705154730 


2.04948272662971 


10" 


-4 


5.25008778789724 


6.30539985918368 


10" 


-2 




6.28719528119678 


7.30407808674598 


10" 


-4 


5.37616000377536 


7.35575817722271 


10" 


-4 




6.38387404484564 


1.68208400736076 


10" 


-5 


5.46963728071208 


1.83331237734851 


10" 


-1 


* 


6.56406612498208 


1.69023199805961 


10" 


-2 


5.48768019361761 


1.77004085558857 


10" 


-6 




6.76964330648490 


1.81072684633499 


10" 


-7 


5.70340946635026 


3.36198337734769 


10" 


-6 




6.79495897876859 


8.47737682200300 


10" 


-3 


5.71149186560478 


3.44485633165334 


10" 


-2 




6.81533915489440 


3.05825161898085 


10" 


-5 


5.78088726970425 


1.07071079267298 


10" 


-4 




6.83002686340033 


9.95790584921162 


10" 


-4 


5.89573570449384 


1.38033053722550 


10" 


-1 


* 



S-(t = 0,p = tt) = 2.52360427892349 S-(t = 0,p = p-) = 5.01384876970501 ■ IP' 1 



LO n (p = TT) 


w n (p = tt) 






Wn(p-) 


w n (p-) 








0.24490318120408 


7.69543336339520 


10" 


-i 


0.87610327625376 


1.95498761012222 


10" 


1 




2.00062423661791 


9.81468201828177 


10" 


-2 


2.50939624323656 


5.59400064382891 


10" 


2 




3.16271478820483 


2.64572507526320 


10" 


-4 


3.47398546669344 


6.94594294796492 


10" 


3 




3.57865017173775 


6.85304507238561 


10" 


-3 


3.60343461539355 


9.85817067330708 


10" 


4 




3.98061972066478 


4.71390118529535 


10" 


-2 


3.71327574076944 


8.01387979043592 


10" 


2 




4.35269425889094 


9.62693115241292 


10" 


-5 


4.17072101766960 


7.78964097133079 


10" 


2 




4.72994233277789 


8.68891912517511 


10" 


-4 


4.31123417306011 


2.56319236222587 


10" 


2 




5.11352567119504 


5.06260574809040 


10" 


-4 


4.77730713351960 


8.29613214115695 


10" 


3 




5.45159687760725 


7.28553124003179 


10" 


-3 


5.12926535498599 


8.40410870562240 


10" 


2 




5.74161870216759 


3.48662020479076 


10" 


-2 


5.41837248921489 


1.83487171330598 


10" 


1 


* 


6.10575246933884 


3.84080214753984 


10" 


-4 


5.66879948906928 


8.77367998473676 


10" 


2 




6.51195982420624 


1.06776716292014 


10" 


-2 


5.94424472415545 


1.41905770467591 


10" 


1 


* 
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TABLE II. The critical quantities 20 (M), r?(M) and 2[1 - a(M)] at M = 1/4 and at the softmode momenta 
p = p- A (M = 1/4) = tt/2, p = p+(M = 1/4) andp = p^{M = 1/4). 



(a) 


Z03[1V1 ) 


V3 \M ) 


2[1 — a:i(p — n/Z, M J J 




p = p 3 (M) 


1.5312 


1.51 


1.54 




(b) 


20i (M) 


??i(M) 


2[l-Q+(p = 7T,M)] 


2[1 - a_ (p = 7r,M)] 


P = 7T 


0.6531 


0.65 


0.62 


0.68 


(c) 


20+ (M) 


<(M) 


2[1 -a+(p = pt(M),M)] 




P = P+{M) 


2.1843 


2.17 


2.40 




(d) 


26»-f (M) 


r?r(M) 


2[l -a-(p = p^(M),M)] 




p = prw 


0.8781 


0.8 - 1.2 


2.1 
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